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We present an analysis of the two-loop amplitude of the muon decay 
in the Standard Model (SM) using algebraic renormalization techniques. 
In addition, we discuss a manifestly BRST invariant IR regulator for the 
photon within the SM. 
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1 Introduction 



In perturbative multi-loop calculations, the subtraction of UV-divergences in quantum 
field theory generally leads to Green functions which fail to respect the symmetries 
of the theory. With the exception of the well-known 75 problem, the method of 
dimensional regularization is compatible with the gauge symmetry but it breaks the 
super symmetry. In these practical method is needed to restore the symmetry 

identities of the gauge symmetry or of the supersymmetry. Here the method of 
algebraic renormalization |lj supplies a complete solution. However, this method has 
rarely been used in practical calculations although it has been applied intensively in 
order to demonstrate the renormalizability of various models. 

In a recent papers 0, 0] , we reviewed the method of algebraic renormalization from 
a practical point of view and proposed an algebraic method combining the advantages 
of the background field method (BFM) and the simplification of (intermediate) Taylor 
subtractions. The method is independent of the regularization scheme; since the local 
breaking terms are under control, one can use the most convenient regularization 
scheme in a specific application. After a straightforward analysis of the corresponding 
(Ward-Takahashi Identities) WTIs and (Slavnov- Taylor Identities) STIs, the spurious 
anomalies introduced by a non-invariant regularization scheme were shown to reduce 
to a few universal breaking terms which depend only on finite Green's functions. The 
method was already applied to several phenomenologically relevant examples in the 
SM, such as the two-loop contributions to the processes H — > 77, to B — > Xgj and 
to the three-gauge boson vertices [0, |, [J. 

Because of the experimental precision of standard model observables at LEP 
(CERN, Geneva), at SLC (SLAC, Stanford) and at TEVATRON (FERMILAB), cal- 
culations of quantum corrections on the two-loop level are necessary; and the 75 play 
a critical role here. The purpose of this note is to offer a theoretical analysis of the 
electroweak two-loop contribution to the muon decay using our algebraic method, 
and to show its efficiency 

Since a detailed self-contained discussion of the fundamental symmetry constraints 
for the SM, of the algebraic renormalization procedure in the BFM, and of our specific 
subtraction method can be found in [Q, [|, we restrict ourselves here to reviewing 
briefly the basic steps of our method (Sec. 2). Then we discuss the muon-decay 
amplitude, in particular the gauge-invariant subset of two-loop diagrams which is 
sensitive to the 75 problem (Sec. 3). It is well-known that there is a physical infrared 
divergence present in the muon-decay amplitude. For this purpose, we propose an IR 
regulator that is manifestly compatible with all the symmetries of the SM (Sec. |J). 
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2 General Strategy 



In the following, we briefly review the main steps elaborated in []2], |3| to renormalize 
a gauge model with a non-invariant regularization technique. The BFM turns out 
to be very important for our purposes and, therefore, we quantize the SM in the 't 
Hooft background gauge || |6|, [/J. 

The use of a non-invariant regularization scheme induces breaking terms into the 
STIs 

S (r (n) ) = rA (n) ' 5 + 0{h n+1 ) , (1) 



which implement the Becchi-Rouet-Stora-Tyutin (BRST) symmetry, and into the 
WTIs 

W (A) (r (n) ) = h n A (n) > w + 0(h n+1 ) , (2) 

which implement the background gauge invariance of the SM. The definition of S and 
W( A ) is given in [§. 

The local breaking terms are denoted by A^' 5, and A^' w . Note that the locality 
is a consequence of the Quantum Action Principle. Here and in the following, de- 
notes the n-loop order regularized and (minimally) subtracted one-particle-irreducible 
(1PI) function. r( n ) includes the renormalization of all subdivergences. The STIs and 
the WTIs are not able to fix the Green functions completely. Indeed it is possible 
to add invariant local terms to the action changing the normalization conditions of 
the Green functions. A complete analysis on the normalization conditions for the SM 
can, for instance, be found in j7|. 

Acting on the broken WTIs (0) with the Taylor operator (1 — T s ) one gets 

(1 - T 6 )W W (r<»>) =0, (3) 



where 5 has to be chosen in such a way that (1 — T 5 )A^' S ^ W = 0. After commuting 

; (A); 



the Taylor operator (1 — T 5 ) with Wm, we obtain 



(A) 



T 5 ')T {n) ] = \t 5 w x) - w (A) T 5 'l r (n) = h n ¥ n) > w (\), (4) 



where 5' is the naive power-counting degree of T^ n \ In general, one has 5 > 5', hence 
the commutation of the Taylor operator with Wm leads to over-subtractions of 
and, thus, to the new breaking terms vJ/W'^A) occurring on the r.h.s. of Eq. @ (for 
more details see 0, |3j). The breaking terms ^^'^(A) for the STIs are defined in the 
same way. Therefore, the application of the Taylor subtraction on Eqs. ([!]) and (0) 
transforms them into 

S (f< n >) = h n ^' s + 0{h n+l ) and w (A) (f (n) ) = n n m^ s + o{h n+1 ) , (5) 
where f» = (1-T 5 )F^. 
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The breaking terms vE^"-)' 5 and ^^' w can be expressed in terms of a linear com- 
bination of ultra-violet (UV) finite Green functions. Here we assumed that up to 
the (n — l)-loop order the Green functions are already correctly renormalized. The 
main difference between vl/M' 5 and ^^>' w is due to the linearity of the corresponding 
operators S and Wm- In the former case one has to consider non-linear terms aris- 
ing from lower loop orders. On the contrary, the linearity of the WTIs enormously 
simplifies the evaluation of the breaking terms and of counterterms. 

Finally, we introduce 

]p(n) _ p(n) _ w(n) _ n _ J^-p(n) _ ^(n) ^ ^ 

where S( n ) is chosen in such a way that the following identities are fulfilled: 

s (r (n) ) = o , W(A) (r (n) ) = o . (7) 

In general, it is quite simple to compute the counterterm, r aT ' = T S T^ + 5^, 
as it can be expressed in terms of Green functions expanded around zero external 
momenta. 

As already mentioned above, there is still the freedom to add invariant counter- 
terms. In other words, we have the freedom to impose normalization conditions that 
lead in addition to Eqs. (|7|) to the following equation being fulfilled: 

M(r {n) )=o, (8) 

where the index i runs over all independent parameters of the SM. As the Green 
function also has to fulfill this condition we have for the counterterm 

Mi (T s r^ + H (n) ) = 0, (9) 

which is a local equation. This means that, whenever the effort to impose the nor- 
malization conditions is done the changes due to the subtraction are only a local 
changes which can be easily compensated. For clarity let us consider an example: for 
the condition on the W boson mass we could choose Mi(T^ + ^_) = T^J^ / _(p*) = 
where the superscript T stands for the transverse part and Re(p*) = Mjy. Notice 
that the imposition of normalization conditions is a very important ingredient for the 
computation in order to compare with other schemes and in order to simplify the 
breaking terms themselves. 

The procedure described so far is heavily based on the Taylor operator T s . In the 
presence of massless particles this may introduce infra-red (IR) divergences. In H we 
presented a modified procedure which resolves this spurious IR problem generally. 
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Figure 1: Example of 0(Nfa 2 ) contributions to the muon decay amplitude with a two- loop 
three-point, two-loop two-point function and with box contribution 

3 Muon Decay Amplitude 

3.1 General settings 

• Muon decay amplitude 

We want to focus on the 0(Nf<y 2 ) contributions including the two-loop three- 
point functions ^w+ Ve (Pu,Pe) Qwith a W and an electron (muon) and electron- 
(muon-) neutrino. This subgroup of 0(Nfa 2 ) contributions to the muon decay 
amplitude is the most delicate one regarding the 75-problem. An example is 
given in Fig. [l]. 

There are further contributions at the two-loop level such as the diagrams in- 
cluding the two- loop-two-point function T^+^-(p) (Fig. [I]). These have already 
been discussed within our approach in Q, but there are no problems with 75 
there. Moreover, there are two-loop box diagrams with a gauge-boson self- 
energy inside (as shown in the last diagram of Fig. [I]). 

Thus, let us focus on contributions like the one shown in the first picture of 
Fig. HJ. We have to consider the complete gauge-invariant subset of two-loop 
contributions to the three-point function Y^+ Ve {p V) p e ). Actually, there are 
various types of diagrams shown in Fig. [J. Here we note that only the first two 
diagrams in Fig. ^| change when switching from a conventional gauge to the 't 
Hooft-background gauge. 



* All momenta are considered as incoming. In the Green functions T<j >1 ...<j >ri they are assigned to 
the corresponding fields starting from the right. The momentum of the most left field is determined 
via momentum conservation. W denotes the background field corresponding to the quantum field 
W. 
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Figure 2: Gauge- invariant subset of 0(Nfa 2 ) contributions to the muon decay including 
vertex and external leg corrections. 



IR problems 

Among the 0(Nf<y 2 ) gauge-invariant subset of diagrams (see Fig. (§)), one has 
to consider those with a virtual photon. Those diagrams are potentially IR 
divergent and only the 0(Nfa 2 ) contributions to the physical amplitude, after 
the inclusion of the Bremsstrahlung radiation of the external electron and muon, 
turn out to be finite. 

Historically ||, the virtual photon contributions, namely the pure QED cor- 
rections, are analysed separately from the electroweak corrections by using a 
suitable decomposition of massless propagators. In that spirit, the pure QED 
corrections have been computed in || |H], [TT[] and the remaining complete elec- 



troweak corrections are IR finite. In paper |L2| the 0(Nf<y 2 ) corrections are 
computed in the MS scheme and the massless quark approximation has been 



used. In [13], [14] the exact fermionic contributions, in on-shell scheme, are taken 
into account. However, since we are not interested in the explicit evaluation of 
the muon amplitude, but only to present a procedure to handle the 75 prob- 
lem in the present process, we will not disentangle the QED corrections in our 
considerations. 

As a consequence, we have to keep the possible IR divergences under control, 
namely we have to be sure that all the steps of the computation are IR regu- 
lated. The situation is worsened by the fact that, according to our procedure, 
the Taylor expansion in Eqs. (^) is performed at zero momentum. For those 
purposes, we propose a BRST invariant IR regularization for the photon within 
the SM (see Sec. f|). This method regulates both physical and spurious IR 
divergences simultaneously. 



5 



An alternative approach to IR problems is the following: regarding 75, the 
delicate diagram is shown in the first picture of Fig. (|1]) which belong to the 
pure vertex corrections T^+ Ve (p v , p e ). Luckily, this vertex Green function is 
the simplest one - compared to box and external-leg corrections - from the IR 
perspective. This is because its physical IR singularities are induced only by 
the on-shell wave-function renormalization. 

Therefore, we can also avoid the physical IR problems by choosing an off-shell 
renormalization procedure. There are two ways of doing this: i) either one 
imposes an on-shell renormalization for the neutrino and, as a consequence of 
WTIs, an off-shell renormalization prescription of the electron is automatically 
provided (see next section), ii) or one can also choose a MS wave-function 
renormalization for external fermions which is infrared finite and compatible 
with the background gauge invariance. 

Finally, to handle the spurious IR problems generated by means of the Tay- 
lor subtraction, a modification of the procedure is discussed in [Q. Here, the 
modified breaking terms ^)W> W 1 occurring in the WTI for the Green function 
Fw+ Ve {Pv,Pe), is written explicitly. 

• Kinematic approximations 

The specific kinematic situation allows for some simplifications: in r^+_ e (py, p e ) 
the W is off-shell, while both the electron and the electron-neutrino are on-shell. 
We can make the approximation pyy = and neglect m muon /mvi/-terms, because 
the muon is almost at rest. All momenta squared are zero, and hence 

PW=Pu=Pe = 0. (10) 

In the following we will derive all symmetry constraints for the general kinematic 
case and then specify to the zero- momentum setting ([L0"D . 

• Subdivergences 

Considering only the 0(Nf<y 2 ) contributions to the two-loop three-point ver- 

(2) 

tex function T ] ^ + _ e (p u ,p e ), we have to take into account three kinds of one- 
loop subdivergences: the three-gauge boson vertices with one background and 
two quantum fields (they have been largely discussed in Q within the BFM 
framework), the quantum gauge boson self-energies (they have been analysed 
in [J without and with the BFM; in particular, in J3] the conversion from 
background field amplitudes to quantum ones is completely exploited) and the 
one-loop T^ + _ (p u ,p e ) amplitude together with their corresponding scalar ver- 
tices where the gauge boson is replaced by the Goldstone boson (notice that 
this amplitude appears also as subdivergence for two-loop three-gauge boson 
vertices and it is extensively discussed in 0]). 
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Moreover, the renormalization of one-loop amplitudes can be quite easily han- 
dled within different regularization techniques, therefore we can assume, for 
the time being, that the one-loop Green functions already satisfy the WTIs (or 
the STIs) and fulfill certain normalization conditions. Nevertheless, to apply 
our procedure, we have to compute the one-loop counterterms which must be 
inserted in one-loop graphs. By using the notation of the introduction, these 
counterterms are given by 

jp(ra) _ p(n) _ ^(n) _ p(n) _ |"j><5p(n) _|_ ^(n) 



bare 



UV 



TT 



bare 



rp5 r (n) , w ( n ) 
1 1 UV ' " 



In the second line we have introduced the bare Green function rl"L in addition. 



bare 



This quantity is defined by r^ n ^ = r b ™j. c — Tyy where Iyy denotes the necessary 
UV counterterms computed in the specific regularization used in the calculation. 
Of course, the complete one-loop counterterms, namely IT*™) — r^ rc , have to be 
taken into account at the two-loop level. 



For instance, in the case of charged gauge-boson self-energies, given T 



(i) 

w+w. 



■ip) 



(which satisfy the normalization conditions and the corresponding STIs) and 



given (p), computed in the same regularization as will be used in the 



w?w, 



two-loop computation of T^ + _ e (p u ,p e ), the counterterms are 



,(1),C.T. 

w+w- 



(p) = Tl r 



,(1) 



■ip) 



(p), 



■ip) 



(12) 



where 



w,i 



W,2 




cfl2 ip(l) 



-<9 2 F (1) 



ip) 
ip) 



— 2 d n (j-d r> i / IT 



a) 



p=0 



p=0 




w+w, 



■ip) 



p=0 



(13) 



In the same way, all the other possible one-loop divergences can be computed 
once the renormalized Green functions are known. 

In some cases, the BFM does not achieve great advantages (for instance, in the 
cases of amplitudes with external fermions only) at the practical level and, on 
the other hand, it could be convenient to use the conventional gauge fixing. 
However, the Green functions computed with external background fields can 
be easily related to those with external quantum fields by using the extended 
versions of the BRST symmetry. 
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3.2 Two-loop vertex function 



Working in the framework of the BFM, there is only one WTI for the vertex function 



W+Ve 



{PviPe) that has to be evaluated at two loops (cf. 



i iPu + Pe) p r^ + _ e (p u , p e ) + i M w Tf + _ e {j> v1 Pe 



+ 



le 



sw 



V2 



T^(-p u )P L -P R Tg>(p 



(2), 



(14) 



Here (p UJ p e ) is a polynomial of the external momenta p v and p e of max- 

imum degree 1. We define the weak mixing angle through the on-shell relation 
Cw = Mw/Mz as we want to maintain the form of the WTIs to be the same to 
all orders. Pl/r — (1 T 7s)/2 are the chiral projectors. 

The breaking terms in ( ]T4|) are generated by a non-invariant regularization pro- 
cedure, for instance by using the 't Hooft-Veltman definition of 75 To remove 
them, according to the procedure described in |2|, [|, we apply the Taylor operator 
(! - T pl, Pe )> obtaining 



i (Pu + Pe 



i 1 T p v , P e) r VK + pI7p e ^' 



Pe 



+ iM } 



t° )r (2) (v v 



.-^-f[fi-TMr^ 



Vu( Pv , 



Pi - Pn 



,Pe 



(15) 



where 



(16) 



are finite and are generated by means of the over-subtraction. 

Notice that the computation of can be also performed without encountering 

any UV divergences. For that purpose, it is sufficient to implement the subtraction 
of UV sub divergences directly in such a way that all the diagrams are always finite. 
Technically, we suggest the use of Zimmermann's subtraction formula ||16|| . Since all 
the integrals involved can be performed analytically, ^x+Ve can ^ e evaluated exactly. 
In this way, we use the Dirac algebra in 4 dimensions without any ambiguities. 

Notice that the zero momentum subtraction in Eq. ( |T5] ) removes exactly the con- 
tribution that we would like to evaluate to compute the muon decay amplitude in the 
approximation (|T0D. However, this contribution can be computed as a counterterm. 
Notice in fact that, the local part of the muon decay amplitude is totally fixed when 
the normalization conditions for fermion two-point functions and the WTIs are used. 
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By using the parametrization 



z/, e 



-^ e ^ e) = C^a^+qW e 7^, (17) 

for two- and three-point functions, and decomposing the breaking terms into scalar 
functions 

we have the solution 



e 

for the coefficients. Notice that the equality ip% — ip\ follows from the consistency 
conditions (see, for example, the discussion in ||) and it provides a check of the 
computation of the breaking terms. In addition, from Eqs. fll9"D , it emerges that 
(f] = i n the case of invariant regularization techniques, namely when tp^ = 
0, Vz. This means that we are not allowed to impose any arbitrary normalization 
conditions for fermion residues. If the neutrino is renormalized in such a way that 
its residues is equal to 1, the electron residue will be clearly different from 1. In this 
way, we have a partial on-shell scheme, we maintain the background gauge symmetry 
and we can avoid the physical IR divergences for the vertex amplitude. 

The final result, namely diagram computation plus counterterms, can be written 
in the following way 



T° F (2) (n nU £ (2) V 1 Pr 4- £ (2) ^ Po 



• (20) 



The parameters ffl and ^ 2 - ) are used to impose mass renormalization conditions on 
the fermion self-energies. 

(2) 

From this we learn that the symmetric amplitude ITV + _ e at zero momentum is 
just given by the universal counterterm: 

r S+,> = ^ = ) = -CUa^-4W M ^" ( 21 ) 

The proposed procedure allows for an efficient computation of the amplitude r?* + _ 

at two-loop order avoiding the 75 problem. In the literature, different techniques 

with different prescription of 75 have been used to evaluate the Feynman diagrams of 

(2) 

IT,j; + _ , however we believe that a rigorous check of these result is desirable. 
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4 Massive U(l) BRST symmetry within the SM 



It is a well-known problem that in the computation of Green's function there are IR 
divergences due to the vanishing photon mass. In this brief section, we will describe 
how to perform a regularization of the photonic IR divergences in a way that is 
consistent with the BRST symmetry of the SM and, thus, preserves the unitarity of 
the model (see also fn|). The choice of such a regulator is motivated essentially by 
the fact that it regulates both the physical and spurious IR divergences. 

We will refer to the Stueckelberg method (see 0, |2(J and references therein). This 
method gives rise to unsolvable renormalization problems in the case of Yang-Mills 
fields, which only can only be resolved through the Higgs mechanism. However, in 
the abelian case it provides a manifestly BRST (and the background gauge) invariant 
model of massive QED. 

For pedagogical purposes, we present a short digression regarding the Stueckelberg 
formalism in QED. The Lagrangian is given by 

£ = --T^F^J 2 + -r-yiAf, d^f + s(aF) + £ ma tter , 

Ag z 2g z m 

^matter = $ ~ + - -dtf)* , (22) 

m 

where F^ u = d^A u - d u A^. ^ 

The apparent non-renormalizable derivative coupling \I/7 M <9 M <£>\I/ can be absorbed 
by a field redefinition 

i> = e^tf, C matter = ^(V^ + Afj) - M)ip. (23) 

Notice that in the non-abelian case a redefinition of fields of type (|23|) will generate 
new non-eliminable non-renormalizable terms p0| . 

The BRST transformations of the fundamental fields are given by 

sA^ = , An =An- —dtf , 

771 

Sip =mc, sc =0, (24) 
sc = b , s b = . 

Thus, the BRST multiplets consist of two trivial pairs (ip,c), (c,b) and one singlet 
Afj, = An — ^dnip- This means that the physical spectrum will be independent of 
(ip, c) and (c, b). 

In the gauge-fixing-ghost term s(cJF), T is a real bosonic function of all the fields 
and their derivatives. For example, the 't Hooft-Feynman gauge fixing 

F=\{b-PA lt -mp), (25) 
g 2 2 
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leads directly to noninteracting ghosts c and c. Moreover, the A — ip sector does 
not contain higher derivatives or dipoles. The field strengths F^ u calculated from A 
and from A coincide. The field b is auxiliary with the algebraic equation of motion 
b = d^A^ + rmp Thus, the C defines a massive abelian gauge field coupled to matter 
fields. 

To extend the Stueckelberg formalism to the SM quantized in the background 
gauge, the scalar field tp transforms under the BRST symmetry^] and under the back- 
ground gauge transformations in the following way 

sip = pic, 5( A )^ = /xA, (26) 

where A is the infinitesimal parameter of the background gauge transformations, /i is 
the IR regulator and c is the U(l) ghost. The latter can be written in terms of the 
combination c = cwCa + s^cz, where ca and cz are the photon ghost and the ghost 
associated with the Z boson, respectively. It follows that a term like 

r stu = J d 4 x f^ffV - nd^B* + y b " b iA > ( 2T ) 

is BRST and background gauge invariant for all the values of the parameter /i and 
the last term provides a mass term for the fields. Other ^-dependent invariant 
terms can be constructed, however it is easy to show that all of them, but (0), can 
be reabsorbed by a simple redefinition of the field B^ (cf. pOl). In addition, in order 
to deal with diagonal two-point functions the 't Hooft gauge fixing 

r stu ' g - f - = s Jd*xc (d^Bn - pn V + £ &) 

= J d 4 x [b (d^B^ -pfnp + £ b) -cd 2 c + pfx 2 c c] , (28) 

is used. Here £o is the conventional gauge fixing parameter of the U(l) sector and p 
is the 't Hooft parameter. With this gauge fixing, it is easy to see that the gauge field 
B^ the scalar field ip and the ghosts c, c (with masses /i 2 , p 2 /i 2 /£o and pfj, 2 ) form a 
quartet which ensures the unitarity of the model. Notice that the BRST variation of <p 
says that this field corresponds to a would-be-Goldstone boson, and the spontaneous 
symmetry breaking mechanism - in the abelian case - can be implemented without 
the Higgs counterpart. 

In the SM framework, the field B^ does not coincide with the physical photon 
field, but the mixing with the third component of SU (2) gauge boson triplet has 
to be considered. We have to cancel this term by modifying the gauge fixing function 
Tb (cf. [0, Eq. (A. 3)) for the abelian field in the following way 



T B = d^B^ + Po ($ + v)H% ($ + v) j + | b — > T B - PW , (29) 



T Notation, definitions, and quantum numbers can be found in 
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where p is the t'Hooft parameter for the p field. By eliminating the Lagrange multi- 
plier b, we have the gauge fixing terms 

= -^[PBr-pw + pofi + v^iQ + v) 

/ 2/2 2 

P V 9 /nun / tt A tt /n\ P g V 



)( ( 9 »B,- m )(HG-HG)-^—- (HG-HG) , (30) 

where g' is the U(l) gauge coupling, v is the vacuum expectation value, G and H 
are the Goldstone boson and the Higgs field, respectively, while G and H are their 
background partners. The first line contains the contribution to the quadratic part 
of the action, this shows that also the masses of the Goldstone boson G are modified 
by the introduction of the Stueckelberg field (p. The mixed terms pd^B^ and Gd^B^ 
are cancelled (in the restricted 't Hooft gauge) by the mixing terms coming from the 
covariant derivatives of the kinetic terms (i.e. from Eq. (|27|) ). Finally the last terms 
describe the interactions between p and the other fields. As can be noticed all the 
interaction terms depend on the background fields. Therefore, the Stueckelberg p field 
can be generated only if the fields G or H appear as external vertices of the amplitude 
or due to the mixing with the neutral Goldstone boson. This is the only difference 
between the Stueckelberg formalism and its application to the SM quantized in the 
background gauge. 
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